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Abstract 

The squark mass-matrix from the soft supersymmetry (SUSY) breaking sector contains a rich 
flavor-mixing structure that allows 0(1) mixings among top- and charm-squarks while being 
consistent with all the existing theoretical and experimental bounds. We formulate a minimal 
flavor-changing-neutral current scheme in which the squark mixings arise from the non-diagonal 
scalar trilincar interactions. This feature can be realized in a class of new models with a horizontal 
U{\)h symmetry which generates realistic quark-mass matrices and provides a solution to the 
SUSY /i-problem. Finally, without using the mass-insertion approximation, we analyze SUSY 
radiative corrections to the H^bc and h°tc couplings, and show that these couplings can reveal 
exciting new discovery channels for the Higgs boson signals at the Tevatron and the LHC. 



PACS numbers: 12.60.-i, 12.15.-y, 11.15.Ex 



1. Introduction 

Weak-scale super symmetry (SUSY) Q, as a leading candidate for new physics beyond the stan- 
dard model (SM), sensibly explains electroweak symmetry breaking, but leaves the understanding 
of flavor sector as a major challenge. Being a new fundamental space-time symmetry between 
fermions and bosons, SUSY necessarily extends the SM flavor structure to include superpartners 
for all fermions and thus adds further puzzles to the flavor physics. To be consistent with the 
experimental data, supersymmetry has to be broken. The breaking of the SUSY lifts the mass 
spectrum of the superpartners and is parametrized by the soft-breaking Lagrangian of the Minimal 
Supersymmetric SM (MSSM) with a large number of free parameters. The soft breaking sector is 
often problematic with low-energy flavor changing neutral current (FCNC) data without making 
specific assumptions about its free parameters. One of the most popular assumptions is the pro- 
portionality of the scalar trilinear A-terms to the fermion Yukawa couplings. This is however not 
a generic feature from the string-theory constructions, and certain forms of non-diagonal A-terms 
and their interesting implications were studied recently J2|, Q . 

In this work, we focus on the flavor-mixings of three-family squarks which originate from the 
scalar mass and the trilinear interaction terms. We observe that the current data mainly suppress 
the FCNCs associated with the first two family squarks (and in some cases, with the first and third 
families), but allow the flavor-mixings between the second- and third- family squarks, the scharm 
(c) and stop (t), to be as large as 0(1) Furthermore, the 0(1) t — c mixings arising from the 
non-diagonal A-term are consistent with the theoretical bounds derived from avoiding the charge- 
color breaking (CCB) as well as maintaining the vacuum stability (VS) ||. Taking a bottom-up 
approach, we first formulate a minimal SUSY FCNC scenario, called the Type-A models, in which 
all the observable FCNC effects come from the non-diagonal trilinear A-term in the c — t sector. 
Then, using the simplest horizontal U(1)h symmetry (d la Froggatt-Nielsen ||), we construct a class 
of new models, called the Type-B models which not only exhibit similar flavor-mixings in the t — c 
sector but also generate realistic quark- mass/mixing pattern and provide a solution to the SUSY 
/i-problem. Such minimal FCNC schemes can reduce the general 6x6 squark-mass-matrix down 
toa4x4or3x3 matrix involving only the c — t sector, and therefore simplify the exact squark 
mass-diagonalization and rotations without invoking the crude mass-insertion approximation. This 
allows quantitative understanding and predictions of the relevant FCNC signatures from the squark 
sector, and thus provides reliable probes to the fundamental SUSY flavor structure encoded in the 
soft-breaking Lagrangian. 

As we will show, exploring the SUSY flavor sector is also important for revealing exciting new 
discovery signatures from the weak-scale supersymmetry, in addition to probing the mechanism 
of soft SUSY breaking. Applying the minimal FCNC schemes, we analyze the SUSY radiative 
corrections to the H^bc coupling and show that this flavor-mixing coupling can be significant to 
provide new discovery signals of charged Higgs via charm-bottom fusion j?]] at the on-going Fermilab 
Tevatron Collider and the CERN Large Hadron Collider (LHC). Then, we further study the flavor 
changing top-decays into the charm-quark and lightest neutral Higgs boson in our schemes, and 
their observability at the LHC. 

2. Minimal Supersymmetric FCNC Models 

2.1. Type-A Minimal SUSY FCNC Models with Non-Diagonal A-Term 

The MSSM soft-breaking squark-sector contains the following quadratic mass-terms and trilinear 
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j4-terms: 

-QliM^Qj - U^M^ijUj - Dl(Ml)ijDj + {A i iQ i H u U j - A^QiH d Dj + c.c.) 
with i,j = 1, 2, 3 being the family indices. This gives a generic 6x6 mass matrix, 
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with m WjZ the masses of (W^, Z°) and M u the up-quark mass matrix. For convenience, we will 
choose the super Cabibbo-Kobayashi-Maskawa (CKM) basis for squarks so that in Eq. (Q), A u 
is replaced by A' u = KuLA u Kjj R and M u by M„ iag , etc, with Kul,r the rotation matrices for 
diagonalizing M u to Mu iag . In our minimal Type- A scheme, we consider all large FCNCs to solely 
come from non-diagonal A' u in the up-sector, and those in the down-sector to be negligible, i.e., we 
define, at the weak scale, 

'0 0' 

A' u = x A, (4) 
I y 1 J 

where, x and y can be of O(l), representing a naturally large flavor-mixing in the t — c sector. 
Such a minimal scheme of SUSY FCNC is compelling as it is fully consistent with the stringent 
CCB/VS theory bounds Q as well as the existing data (In the case that the CKM matrix is 
generated from the down-quark sector alone Q, A' u simply reduces back to A u .) Similar pattern 
may be also defined for A^ in the down-sector, but the strong CCB/VS bounds permit 0(1) b — s 
mixings only for very large tan/3 because m& <C mt- Thus, to allow a full range of tan/3 value, we 
consider an almost diagonal A^. Moreover, indentifying the non-diagonal A u as the only source of 
the observable FCNC phenomena in the Type-A models implies the squark-mass-matrices M~ ~ in 
Eqs. (||)-(||) to be nearly diagonal. For simplicity, we define 
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with fho a common scale of scalar-masses Q . 

Within this minimal Type-A scheme, we observe that the first family squarks ul,r decouple 
from the rest in (||) so that the 6x6 mass-matrix is reduced to 4 x 4, 
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In Eq. (H), we ignore tiny terms of 0(m c ) or smaller. The reduced squark mass matrix (JgJ) has 6 
zero-entries in total and is simple enough to allow an exact diagonalization. Especially, for the cases 
of (i) x ^ 0, y = 0, (called Type-Al) and (ii) x = 0, y ^ 0, (called Type-A2), one more squark, cr 
(in Type-Al) or cl (in Type-A2), decouples, and the Eq. (||) further reduces to a 3 x 3 matrix. We 
have worked out the general diagonalization of 4 x 4 matrix (||) for any (x, y). The mass-eigenvalues 
of the eigenstates (ci, c 2 , ti, £2) are, 
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where u± 
sector as 



X? + (A x ± A y ) 2 



From (0), we can deduce the mass-spectrum of the stop-scharm 



M h < M51 < M- c2 < M- t2 . (9) 

In Eq. (||), the stop i\ can be as light as about 100—300 GeV for the typical range of m > 0.5—1 TeV. 
The 4x4 rotation matrix of the diagonalization is given by, 



(10) 



with 









C1C3 C1S3 S1S4 


S1C4 




/ - \ 

Cl 


CR 






C2S3 c 2 c 3 s 2 c 4 


-S2S4 




C2 


h 






S1C3 -S1S3 cis 4 


cic 4 




h 


, tR , 




*. 


•S2S3 -S2C3 c 2 c 4 


-C2S4 ( 








1 

"71 




1/2 

, s 4 = 


1 


or, s 3 = 


= l/v2 (if xy ^ 0), where s 


2 , „2 _ 

i + c i ~ 


1 





(11) 



and S3 = (if xy = 

2.2. Type-B Minimal SUSY FCNC Models with a Horizontal U(l) Symmetry 

The minimal Type-A SUSY FCNC schemes with a non-diagonal A u -term, as discussed above, 
are truly economical as they uniquely result from imposing all the stringent theoretical and exper- 
imental bounds. Below, we further support such type of FCNC in the t — c sector by providing 
theoretically compelling constructions based upon a minimal family symmetry. An attractive ap- 
proach is to make use of the simplest horizontal U(1)h symmetry to generate realistic flavor structure 
of both quarks and squarks, via proper powers of a single suppression factor || [|, and to provide 
a solution to the SUSY //-problem. For convenience, we define the suppression factor e = (S)/A 
to have a similar size as the Wolfenstein-parameter A in the CKM matrix, i.e., e ~ A ~ 0.22 0. 
Here, (5) is the vacuum expectation value of a singlet scalar 5, responsible for spontaneous U{1)h 
breaking, and A is the scale at which the U(1)h breaking is mediated to light fermions. In general, 
the supermultiplets of three- family quarks/squarks may carry different U(1)h charges, as defined 
in Table 1. 



Table 1. Quantum number assignments under the horizontal symmetry U{1)h- 
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It is then straightforward to deduce the following hierarchy structures in the quark mass matrices: 

(12) 
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Similarly, in the CKM matrix, 

(K s , Vk, Kt) ~ (a' 11 '^, X h2 ~ hs , X h ^) 



(13) 



Unlike Ref. ||, the key ingredient of our model-buildings is to impose a new condition, 

a 2 = a 3 , (14) 

which ensures the mixings between c and t in the squark mass matrix to be naturally of O(l). 
From the condition ([lj) and the current data of quark-masses and CKM angles (which can all be 
counted in powers of A), we find an almost unique solution for all quark/squark quantum numbers 
(cf., Table 2), which is based upon a single U{V)h symmetry and will be called the minimal Type-B 
scheme hereafter. In Table 2, we consider tan/3 ~ O(l) for the down-type quark mass-matrix 
[cf. Eq. (O)]. The extension to a larger tan/3 only affects the quantum numbers of d/s in a trivial 
way as it just contributes an overall factor 1/tan /3 ~ X k (with the integer k ~ 0.66 log tan /3) to 
in Eq. (12) and thus simply adds — k to each quantum number of dj listed in Table 2. 



Table 2. Quantum number assignments are derived for the minimal Type-B model with tan (3 ~ 0(1). 
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There are some slight variations of this minimal Type-B model, by allowing the quantum num- 
bers of Qj's to have £(£')-dependence, but they all predict the same patterns for the quark masses 
and mixings. By this construction, we have attempted to simultaneously solve the SUSY ^-problem 
from the same U (1)h- A dynamical //-term can originate from (k/ A n ~ 1 )S n H u H ( [ , with n = £ + 
such that n = K\ n ^ 1 (S) is generated at a scale (5) <C Mpi anc k. Hence, a weak-scale value of /i can 
be obtained by properly choosing n for a given (5). If U(1)h is not related to the SUSY //-problem, 
the minimal Type-B model becomes truly unique, corresponding to a special case of £ = £' = 
in Table 2. With Table 2, we can readily derive the structures of quark/squark mass- matrices. For 
instance, the up-quark mass-matrix takes the form of 

f A 7 A 4 A 4 ^ 
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for any tan/3 > 1, while the squark mass-matrices M| L and M^ fi in Eq. @ are deduced as, 
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Eqs. (|l5|) and ( |i~6| ) show a partial quark-squark "alignment" that effectively suppresses the FCNCs 
between 1st and 2nd(3rd) families, while at the same time provides 0(1) mixings in the t — c sector 
of M RR . Furthermore, because squarks carries the same t/(l)#-charge as quarks, the t — c mixings 
originating from a non-diagonal A u term is predicted to share the same hierarchy structure as that 
in M u [cf. Eq. (|i~5|)1. This, however, does not imply an exact "proportionality" between A u and M u 
because the power-counting of A allows the coefficients in the corresponding entries of A u and M u 
to differ by ratios of 0{1). Ignoring the small 0(A 3 ) ~ 1% terms in M u , we can diagonalize M u 
by a 2 x 2 rotation of the singlet quarks (c, t). Under the super-CKM basis, the off-diagonal block 
M\ R of the Eq. @ becomes 



M 2 LR = A' u v sin p/y/2 - M^ iag pcatp. 



(17) 
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A u Ky R + 0(A 3 ) and the singlet-quark rotation matrix Kjjr contains a 



where A' u 

nontrivial sub- matrix involving only the second and the third family squarks. After neglecting the 
tiny 0(A 3 ) terms, we can parametrize the minimal A' u term of the Type-B model as 
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where the size of the mixing parameter y is naturally of O(l). Thus, under this construction, the 
squarks (ul, ur, c£) decouple and the Eq. (0) greatly reduces to a 3 x 3 matrix, which takes the 
form, under the basis (cr, t^, tn), 
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In the above equation, A y = yAv sin/?/\/2, and the parameter x = O(l) characterizes the mixing 
of cr — Ir in the mass matrix M RR [cf., Eq. fll6|)]. For convenience, we define the typical case with 
y = as the Type-B 1 scheme, and that with x = as the Type-B2 scheme. We note that the 
Type-B2 model is identical to the Type-A2 model as they have the same form of the non-diagonal 
A u and ex decouples. Also, Type-B2 mass matrix with y = in Eq. (|l9|) takes the same structure 
as that of Type-Al except that Eq. (|l~9| ) involves cr (not c£) and its x originates from M RR (not 
Au). Without losing generality, we will study the physical applications of the typical Type-Al and 
-A2 (-B2) models in the following sections. 



3. SUSY Radiative Corrections and Higgs Signatures at Colliders 

3.1. SUSY Induced H^bc Vertex and Production at Hadron Colliders 

Different from the TopColor models and the Type-Ill of two-Higgs-doublet models discussed in 
Ref. Q, the MSSM Higgs sector has no FCNC at the tree level and the same is true for flavor- 
changing mixings in the charged sector (except the usual CKM mechanism). Thus, the new flavor- 
changing effects in the neutral and charged sectors of the MSSM have to be generated radiatively. 
From Eq. ( |l0| ) and its resulted Feynman rules, we can compute the dominant SUSY-QCD radiative 
corrections to the H^bc coupling. It contains vertex corrections [scharm(stop)-sbottom-gluino loop] 
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and self-energy corrections [scharm(stop)-gluino loop], and can be summarized as, 



Y H+h - c = iu c {k 2 ) (F L P L + F R P R ) u 6 (fci) , 
Fr p - F° + F v 4- F s ^ 
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where Pl,r = (1 T 7s)/2 and the tree-level results are 



(F°, F° ) = -^*- (m c cot £ m b tan /?) , (21) 

with V^fe being the CKM mixing matrix element involving c and b quarks. The one-loop vertex 
corrections in the Type-Al model give, 



f r = ^mg^2Kf k C (m 2 H , 0,0; m~ b .,m § ,m % ), 



(22) 



where % G (£2,^1,^2)) &j G {^1^2), and Co is the 3-point C-function of Passarino-Veltman [^]. kQ. 
is the product of the relevant H^-bj-Uk, bj-g-b and u k -g-c couplings. The bi — 6r mixings are also 
included, which can be sizable for large tan (3. Furthermore, the Type-Al self-energy corrections 
yield 



F[ = 0. 



F S = F0^^^(-)i+^ (0;^,m fi , 

* 5=1,2 
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where Bq is the 2-point Passarino-Veltman function [^] and s± is given in Eq. ( |1C|) for y = 0, and 
tree-level H^tb couplings i 7 ^) = (gVtb/ V2m w )(mtcot (3, m&tan/?). In Eqs. (|22|)-(p3|), the tiny 
sub-leading terms suppressed by the powers of m c /m t: g have been ignored. 

The form factors F^' S and F^' S for the Type-A2 and Type-B2 models can be easily obtained 
from Eqs. (p2|)-(|2"3|) by the exchanges of L <-> R and x —> y. We note that Fl (Fr) vanishes in 
Type-Al (Type-A2 & -B2) schemes because cr (cl) decouples. 

We find that the effective H^bc couplings Fl,r are typically around 0.03 — 0.1 for (x, y) ~ 
0.5 — 0.9, (A, mo) — 0.5 — 2 TeV, and tan (3 ~ 15 — 50. Therefore, the production of a charged Higgs 
boson via the b-c quark fusion process can become important in a wide range of the parameter space. 
In Fig. 1(a), the production cross sections of via pp/pp — > H^X at the Tevatron and LHC are 
shown for (//, m s , m ) = (300, 300, 600) GeV, (A, -A b ) = 1.5 TeV, tan/3 = (15, 50), and x = 0.75 
in the Type-Al model. The dotted and dash curves are the SM production rates of cs — > and 
cb — > (induced by V c b — 0.04), respectively, after including the complete next-to-leading order 
(NLO) SM-QCD corrections plf| . (The NLO SM-QCD corrections include the subprocesses with one 
single gluon in the initial state.) The solid curve is the full cb — > production rate as a function 
of the Higgs boson mass after including both the SM-QCD and new SUSY-QCD corrections. As 
indicated, the SUSY loop corrections can significantly dominate over the CKM-suppressed F® R 
contributions by a factor of ~2 — 5. In Fig. 1(b) and (c), we show the X-factor, defined as the ratio 
of (-F| + Fx) over (F® + F R ), which characterizes the enhancement of the production rate by 
the SUSY loop contributions. 
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To study the detection of a scalar, we consider its decay modes in the following. For 
m H ^ 190 GeV, the rv channel dominates, and for mu > 190 GeV, mostly decays into the tb 
channel unless mass is above the threshold of W^hP, in which case the W^hP channel (with 
W — > Iv and h — > bb) can become important as well Fig. 1 suggests that the Tevatron may be 
sensitive to the signals for ran below ~ 300 GeV, with an integrated luminosity of 2 — 20 fb -1 
per detector, while the LHC can potentially probe the full mass-range of with an integrated 
luminosity of about 100 fb -1 . A detailed Monte Carlo simulation is needed to further quantify the 
discovery potential of the Tevatron and the LHC, which is beyond the scope of this Letter. 

We have also examined the similar SUSY-induced enhancement K-factor for the production 
rate in the type-A2 and -B2 models and found that in the low tan (5 (< 5 — 10) region it can reach 
to about 2 — 5 for y = 0.5 — 0.9. But, as tan (3 increases to above ~ 15, the enhancement factor 
decreases to less than ~ 1.3 for Higgs mass below 1 TeV. 

3.2. SUSY Induced h°tc Vertex and Neutral Higgs Signal from Top Decay 

It is known that the SM branching ratio of the flavor-changing top decay t — > ch° is extremely 
small (< 10~ 13 - 10~ 14 |2|), so that this channel becomes an excellent window for probing new 
physics [13, 14, 15|. Our minimal SUSY FCNC models predict the branching ratio of the t — ► chP 



decay to be substantially above the SM value so that this decay mode becomes observable at the 
LHC. (Note that this decay channel is always kinematically allowed in the MSSM.) 

In our minimal FCNC schemes, the one-loop SUSY QCD induced tchP coupling can be written 

as 

Ttch = iu c (k 2 ) (F L P L + F R P R ) u t (h) , 

(24) 

'L.R T r L,R ) 



Fl,r = Fl R + F[ 



which contains the vertex corrections (from scharm-stop-gluino, stop-stop-gluino and scharm-scharm- 
gluino, triangle loops), and the self-energy corrections (from stop-gluino and scharm-gluino loops). 
The one-loop vertex corrections in Type-Al are, 

F L = ~^J2 X f kmt ( C ° + C ii){ m hi m u^ m uji m ~gi m u k ) , 

(25) 

j,k 

where Uj^ G (£2,^1,^2), and (Co,Cn) are the 3-point C-function of Passarino-Veltman. is 
the product of the relevant h — Uj — Uk and Uk — g — t(c) couplings, derived from applying the 
squark-rotation (|l0|). The Type-Al self-energy corrections yield, 

Fl = 0, 

(26) 

Fr = Fo^^[B (0-,m- g ,m i2 )-B (0-rn- g ,m h )] , 



where Bq is the 2-point Passarino-Veltman function and sg is given in Eq. ( |i~0P with y = 0. Fq 
denotes the tree- level h° — t — i coupling, and is given by F° = (mt/v) (cos a/ sin/3). Again, in 
Eqs. (p^)-(p^), we have ignored the tiny sub-leading terms suppressed by the powers of m c /m tt g. 
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The form factors F L ' R in the Type-A2 model can be obtained from (]25|)-(p6[) by the exchanges 
of L <-> R and x — * y everywhere. 

For the numerical study, we assume that the only dominant decay mode of the top quark is its 
SM decay mode, t — ► bW . Thus, the decay branching ratio of t — ► ch° is given by, Br[t — > c/i°] ~ 
T[i -> c/i°]/r[t 6W], with the partial decay width 

i 

2 (Fl + F 2 R ), (27) 

where the form factors (Fl, Fr) are defined in Eq. ([24]). As summarized by Table 3, in our minimal 
SUSY-FCNC schemes, the decay branching ratio Br[i — > c/i°] can be as large as 10 -3 — 10 -5 over 
a large part of the SUSY parameter space where the mass of the lightest Higgs boson h° is around 
110— 130 GeV. We see that these decay branching ratios are very sensitive to the mixing parameter 
x when it varies from 0.5 to 0.9. One reason is that the branching ratio (or decay width) contains, 
besides other mass-diagonalization effects, a sensitive overall power factor x 2 which comes from the 
stop-scharm mixing induced flavor-changing coupling in the squark-squark-gluino triangle loop and 
squark-gluino self-energy loop. Another reason is that unlike the usual analyses with mass-insertion 
approximation, we have performed exact squark mass diagonalization [cf. Eqs. (8)-(ll)], so that 
stops and/or scharms can have significant mass-splittings. For instance, Eq. (8) shows that the 
mass-splitting between two top-squarks is not just due to the usual left-right mixing from X t , but 
also arises from the non-diagonal A-term, xA (and yA). The latter further enhances the mass- 
splittings and results in a light t\ of mass ~ 300 — 100 GeV for x = 0.5 — 0.9, and the heavier t2 
always has mass above the input mo (set as 600 GeV in Table 3). The radiative vertex corrections to 
tch° coupling are thus dominated by the diagram with t\ — t\ — g triangle loop as we have explicitly 
verified from the 3-point Passarino-Veltman C-functions in Eq. (|25|). This second reason further 
increases the sensitivity of our decay branching ratios to the mixing parameter x. From Table 3, 
we also see that for moderate mixings with x < 0.5, the branching ratios are generally bounded to 
around the order of 10 -5 , consistent with other studies in the literature [O]. The similar conclusion 
also holds for our Type-A2 and -B2 models. 

Since the LHC with an integrated luminosity of 100 fb _1 can produce about 10 8 t and t pairs [jOj ], 
it can have a great sensitivity to discover this decay channel and test the model predictions, by 
demanding one top decaying into the usual bW^ mode and another to the FCNC ch° mode. As 
shown in a recent model-independent Monte Carlo analysis p4j , the LHC (100fb _1 ) can already 
measure the Br[t -> ch°] down to the level of 4.5 x 10" 5 at the 95% C.L. The future Linear Collider, 
with a high luminosity, is also expected to have a good sensitivity to detect this channel. 

Table 3. Br[t — ► ch°] x 10 3 is shown for a sample set of Type-Al inputs with (mo, fJ>,A) = 
(0.6, 0.3, 1.5) TeV and Higgs mass M^o = 0.6 TeV. The three numbers in each entry correspond 
to x = (0.5, 0.75, 0.9), respectively. 
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4. Conclusions 



The three-family squark mass-matrix originating from the soft SUSY breaking sector contains 
a rich flavor-mixing structure. In this work, we have constructed the minimal FCNC schemes for 
the squark mass-terms and the scalar trilinear interactions which are consistent with the existing 
experimental and theoretical bounds. We find that the 0(1) large mixings among the top- and 
charm-squarks are allowed. We demonstrate that this feature can be naturally realized in a class 
of new models with a horizontal U{\) symmetry which also generate realistic quark- mass pattern 
and solve the SUSY /^-problem. Finally, we systematically analyze the dominant super symmetric 
radiative corrections to the bcH^ and tchP couplings in our minimal schemes, without using the 
mass-insertion approximation. We show that these couplings can be significant to provide new 
discovery signatures of the charged and neutral Higgs bosons at the Fermilab Tevatron and the 
CERN LHC. 
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Figure 1: (a). production via cb (and cs) fusions at hadron colliders, with sample inputs 

tan/3 = 15(50) shown as lower (upper) set of curves, and x = 0.75. (b) and (c). Factor K = 
(F£ + Fft) I [F® 2 + Fft 2 J for H^bc vertex, as a function of parameter x and for tan/3 = (15, 50). 
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